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Abstract 

Our talk at Lisbon SAMP conference was based mainly on our recent results on small 
diameter asymptotics for solutions of the Helmgoltz equation in networks of thin fibers. These 
results were published in |21] . The present paper contains a detailed review of |21] under 
some assumptions which make the results much more transparent. It also contains several new 
theorems on the structure of the spectrum near the threshold, small diameter asymptotics of 
the resolvent, and solutions of the evolution equation. 
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1 Introduction 

The paper concerns the asymptotic spectral analysis of the wave problems in systems of wave guides 
when the thickness e of the wave guides is vanishing. In the simplest case, the problem is described 
by the stationary wave (Helmholtz) equation 

— e 2 Au — Ait, x € f2 e , Bu = on <9f2 e , (1) 

in a domain fl £ C R , d > 2, with infinitely smooth boundary (for simplicity) which has the following 
structure: fi £ is a union of a finite number of cylinders Cj_ E (which we shall call channels), 1 < j < N, 
of lengths lj with the diameters of cross-sections of order O (e) and domains Ji j6 , • • • , Jm, £ (which 
we shall call junctions) connecting the channels into a network. It is assumed that the junctions have 
diameters of the same order 0(e). The boundary condition has the form: B — 1 (the Dirichlet BC) 
or B = ^ (the Neumann BC) or B = + a(x), where n is the exterior normal and the function 
a > is real valued and does not depend on the longitudinal (parallel to the axis) coordinate on the 
boundary of the channels. One also can impose one type of BC on the lateral boundary of f2 e and 
another BC on free ends (which are not adjacent to a junction) of the channels. For simplicity we 
assume that only Dirichlet or Neumann BC are imposed on the free ends of the channels. Sometimes 
we shall denote the operator B on the lateral surface of f2 e by Bo, and we shall denote the operator 
B on the free ends of the channels by B e . 

Let m channels have infinite length. We start the numeration of C Ji£ with the infinite channels. 
So, lj = oo for 1 < j < m. The axes of the channels form edges Tj, 1 < j < N, of the limiting 
(e — > 0) metric graph T. We split the set V of vertices Vj of the graph in two subsets V = V\ U V 2 , 
where the vertices from the set V\ have degree 1 and the vertices from the set V<i have degree at least 
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Figure 1: An example of a domain Q £ with four junctions, four unbounded channels and four 
bounded channels. 

two, i.e. vertices Vi € V\ of the graph T correspond to the free ends of the channels, and vertices 
Vj G V2 correspond to the junctions Jj e . 

Equation degenerates when e — 0. One could omit e 2 in (Q]). However, the problem under 
consideration would remain singular, since the domain fl £ shrinks to the graph r as e — > 0. The 
presence of this coefficient is convenient, since it makes the spectrum less vulnerable to changes in e. 
As we shall see, in some important cases the spectrum of the problem does not depend on e, and the 
spectrum will be proportional to e~ 2 if e 2 in ([T]) is omitted. The operator in L 2 (fl £ ) corresponding 
to the problem ([!} will be denoted by H £ . 

The goal of this paper is the asymptotic analysis of the spectrum of H £ , the resolvent (H £ — A) -1 , 
and solutions of the corresponding non-stationary problems for the heat and wave equations as s — > 0. 
One can expect that H £ is close (in some sense) to a one dimensional operator on the limiting graph 
r with appropriate gluing conditions at the vertices v € V. The justification of this fact is not always 
simple. The form of the GC in general situation was discovered quite recently in our previous paper 

An important class of domains O e are self-similar domains with only one junction and all the 
channels being infinite. We shall call them spider domains. Thus, if f2 e is a spider domain, then 
there exist a point x = x(e) and an e- independent domain f2 such that 

fl £ = {{x + ex) : x (2) 

Thus, f2 e is the e-contraction of f2 = fij.. 

For the sake of simplicity we shall assume that f2 e is self-similar in a neighborhood of each 
junction. Namely, let Jj( v ).e be the junction which corresponds to a vertex v S V of the limiting 
graph r. Consider a junction Jj( v ) t£ and all the channels adjacent to Jj( v ), e - If some of these channels 
have finite length, we extend them to infinity. We assume that, for each v € V, the resulting domain 
il Vt s which consists of the junction Jj( v ),e an d the semi-infinite channels emanating from it is a spider 
domain. We also assume that all the channels C^ £ have the same cross-section uj £ . This assumption 
is needed only to make the results more transparent. From the self-similarity assumption it follows 
that u> £ is an e— homothety of a bounded domain u> S R d ~ l . 

Let Ao < Ai < A2... be eigenvalues of the negative Laplacian — A^-i in u with the BC Bqu — on 
dui where we put e — 1 in B , and let {p n (y)}, V € w e be the set of corresponding orthonormal 

eigenfunctions. Then A ra are eigenvalues of — e 2 A c ;_i in uj £ and {tfiniy/s)} are the corresponding 
eigenfunctions. In the presence of infinite channels, the spectrum of the operator H £ consists of 
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an absolutely continuous component which coincides with the semi-bounded interval [Ao,oo) and a 
discrete set of eigenvalues. The eigenvalues can be located below Ao and can be embedded into the 
absolutely continuous spectrum. We will call the point A = Ao the threshold since it is the bottom 
of the absolutely continuous spectrum or (and) the first point of accumulation of the eigenvalues 
as e — > 0. Let us consider two simplest examples: the Dirichlet problem in a half infinite cylinder 
and in a bounded cylinder of the length I. In the first case, the spectrum of the negative Dirichlet 
Laplacian in Q £ is pure absolutely continuous and has multiplicity n + 1 on the interval [A„, oo). In 
the second case the spectrum consists of the set of eigenvalues A„. m = A n + e 2 m 2 /I 2 , n > 0, m > 1. 

The wave propagation governed by the operator H £ can be described in terms of the scattering 
solutions and scattering matrices associated to individual junctions of il £ . The scattering solutions 
give information on the absolutely continuous spectrum and the resolvent for energies in the bulk 
of the spectrum (A > A ). The spectrum in a small neighbourhood of Ao and below Ao is associated 
with the parabolic equation. However, the scattering solutions allow us to approximate the operator 
H £ , e — > 0, by a one dimensional operator on the limiting graph for all values of A > Ao- In particular, 
when A w Ao the corresponding GC on the limiting graph are expressed in terms of the limits of the 
scattering matrices as A — > Ao- 

The plan of the paper is as follows. Next section is devoted to historical remarks. A more detailed 
description of the results from our paper [21] on the asymptotic behavior of the scattering solutions 
(A > Ao + S, e — ► 0) is given in section 3. In particular, the GC on the limiting graph are described 
(A > Ao + 8). The Green function of the one dimensional problem on the limiting graph is studied 
in section 4. The resolvent convergence as e — * is established in section 5 when A is near Ao- This 
allows us to derive and rigorously justify the GC for the limiting problem with A w Ao which where 
obtained earlier [5T] only on a formal level, as the limit as A — > Ao of the GC with A > Ao- A detailed 
analysis of these GC is given. 

Let A = Ao + 0(e 2 ). It has been known (see references in the next section) that for an arbitrary 
domain fl £ and the Neumann boundary condition on the GC on the limiting graph L is 

Kirchhoff's condition. The GC are different for other boundary conditions on <9f2 e . It was shown in 
[2T] that for generic domains fi e and the boundary conditions different from the Neumann condition, 
the GC at the vertices of T are Dirichlet conditions. It is shown here that, for arbitrary domain f2 e , 
the GC at each vertex v of the limiting graph has the following form. For any function ? on F, we 
form a vector whose components are restrictions of q to the edges of F adjacent to v. The GC 
at v, with A near Ao, are the Dirichlet condition for some components of the vector <^} v ' and the 
Neumann condition for the remaining components where is a rotation of fW. 

Note that the resolvent convergence provides the convergence of the discrete spectrum. In the 
presence of finite channels Cj iE , the operator H £ has a sequence of eigenvalues which converge to 
Ao as e — > (see the example above). Thus, these eigenvalues are asymptotically (e — > 0) close to 
the eigenvalues of the problem where the junctions are replaced by Dirichlet/Neumann boundary 
conditions. The final result concerns the inverse scattering problem. The GC of the limiting problem 
depend on A if A > Ao + 8. A A-independent effective potential is constructed in the last section of 
the paper which has the same scattering data as the original problem. This allows one to reduce 
the problem in Q £ to a one dimensional problem with A-independent GC. 

2 Historical remarks. 

Certain problems related to the operator H £ have been studied in detail. They concern, directly or 
indirectly, the spectrum near the origin for the operator H £ with the Neumann boundary condition 
on dn £ , see [5JE1 Q21 [13l QUI QU [H [25l [26]. The following couple of features distinguish the Neumann 
boundary condition. First, only in this case the ground states ipo(y/e) = 1 on the cross sections of 
the channels can be extended smoothly onto the junctions (by 1) to provide the ground state for 
the operator H £ in an arbitrary domain Q £ . Another important fact, which is valid only in the case 
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of the Neumann boundary conditions, is that Ao = 0. Note that an eigenvalue A = (i of the operator 
H £ contributes a term of order e - ^ to the solutions of the heat equation in fl £ . The existence of the 
spectrum in a small (of order 0(e 2 )) neighborhood of the origin leads to the existence of a non-trivial 
limit, as e — > 0, for the solutions of the heat equation. Solutions of the heat equation with other 
boundary conditions are vanishing exponentially as e — > 0. 

The GC and the justification of the limiting procedure e — ► when A is near Ao = and the 
Neumann BC is imposed at the boundary of f2 e can be found in [12], [18], [19], [26]. Typically, the 
GC at the vertices of the limiting graph in this case are: the continuity at each vertex v of both the 
field and the flow. These GC are called Kirchhoff 's GC. The paper [T3] provides the convergence, as 
e — > 0, of the Markov process on il £ to the Markov process on the limiting graph for more general 
domains Q £ (the cross section of a channel can vary). In the case when the shrinkage rate of the 
volume of the junctions is lower than the one of the area of the cross-sections of the guides, more 
complex, energy dependent or decoupling, conditions may arise (see [15], [19], [5] for details). 

The operator H £ with the Dirichlet boundary condition on dQ £ was studied in a recent paper 
[24] under conditions that A is near the threshold Ao > and the junctions are more narrow than 
the channels. It is assumed there that the domain f2 e is bounded. Therefore, the spectrum of the 
operator |T]) is discrete. It is proved that the eigenvalues of the operator (Q]) in a small neighborhood 
of Ao behave asymptotically, when e — > 0, as eigenvalues of the problem in the disconnected domain 
that one gets by omitting the junctions, separating the channels in f2 e , and adding the Dirichlet 
conditions on the bottoms of the channels. This result indicates that the waves do not propagate 
through the narrow junctions when A is close to the bottom of the absolutely continuous spectrum. 
A similar result was obtained in P] for the Schrodinger operator with a potential having a deep 
strict minimum on the graph, when the width of the walls shrinks to zero. It will be shown in this 
paper, that the same result (the GC on the limiting graph is the Dirichlet condition if H £ is the 
operator with the Dirichlet boundary condition on dtt £ and A w Ao > 0) is valid for generic domains 
fl £ without assumptions on the size of the junctions. 

The asymptotic analysis of the scattering solutions and the resolvent for operator H £ with 
arbitrary boundary conditions on dfl £ and A in the bulk of the absolutely continuous spectrum 
(A > Ao + S) was given by us in [20], [21]. It was shown there that the GC on the limiting graph 
can be expressed in terms of the scattering matrices defined by junctions of £l £ . Formal extension of 
these conditions to A = Ao leads to the Dirichlet boundary conditions at the vertices of the limiting 
graph for generic domains tt £ . Among other results, we will show here that the asymptotics obtained 
in [50], [H] are valid up to A = Ao- 

There is extended literature on the spectrum of the operator H £ below the threshold Ao (for 
example, see [3]-|ll] and references therein). We shall not discuss this topic in the present paper. 
Important facts on the scattering solutions in networks of thin fibers can be found in [23] , [22] . 

3 Scattering solutions. 

We introduce Euclidean coordinates (t, y) in channels Cj i£ chosen in such a way that i-axis is parallel 
to the axis of the channel, hyperplane -R^ -1 is orthogonal to the axis, and Cj >e has the following 
form in the new coordinates: 

Cj. £ = {(t,ey) :0<t<lj, ye lu}. 

Let us recall the definition of scattering solutions for the problem in fl £ . In this paper, we'll need 
the scattering solutions only in the case of A € (Ao, Ai). Consider the non- homogeneous problem 

(-e 2 A - A)u = /, xen £ ; Bu = on dfl e . (3) 
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Definition 1 Let f £ L 2 com (fL e ) have a compact support, and Xq < X < Ai . A solution u of is 
called outgoing if it has the following asymptotic behavior at infinity in each infinite channel C'j. £ , 
1 < j < m: 

M = a J e l ^ Zi Vo(2//e) + 0(e-^), 7 = 7 ( £ ,A)>0. (4) 
If X < Ao, a solution u of (0) is called outgoing if it decays at infinity. 

Definition 2 Let Ao < X < X\. A function ^ = ^ p e \ < p < m, is called a solution of the 
scattering problem in f2 e if 

(-A - A)* = 0, x e Q e ; = on dfl E , (5) 

and has the following asymptotic behavior in each infinite channel Cj j£ , 1 < j < m : 



*£° = Srje-i^-^Wv/e) + W'^Vofo/e) + 0( e -^), (6) 
where 7 = 7(6, A) > 0, and 8 p j is the Kronecker symbol, i.e. 5 Pl j = 1 if p = j, 5 Pl j =0 if p =/= j. 

The first term in ([6]) corresponds to the incident wave (coming through the channel C Pj£ ), and 
all the other terms describe the transmitted waves. The transmission coefficients t p j = t p j(e,X) 
depend on e and A. The matrix 

T = [tpj] (7) 

is called the scattering matrix. 

The outgoing and scattering solutions are defined similarly when A € (A n ,A n +i). In this case, 
any outgoing solution has n + 1 waves in each channel propagating to infinity with the frequencies 
VA — X s /e, < s < n. There are 771(71+ 1) scattering solutions: the incident wave may come through 
one of m infinite channels with one of (n + 1) possible frequencies. The scattering matrix has the 
size m(n + 1) x 771(77 + 1) in this case. 

Theorem 3 The scattering matrix T, X > Ao, A ^ {Aj}, is unitary and symmetric (t p j = tj tP ). 

The operator H E is non-negative, and therefore the resolvent 

Rx = {H £ -X)- 1 : L 2 (n e ) ^ L 2 {il £ ) (8) 

is analytic in the complex A plane outside the positive semi-axis A > 0. Hence, the operator 
is analytic in k in the half plane Imfc > 0. We are going to consider the analytic extension of the 
operator R k 2 to the real axis and the lower half plane. Such an extension does not exist if R^ is 
considered as an operator in Z/ 2 (fi £ ) since R^ is an unbounded operator when A = k 2 belongs to 
the spectrum of the operator R\. However, one can extend R^ analytically if it is considered as an 
operator in the following spaces (with a smaller domain and a larger range): 

R k i : L 2 com (n £ ) - L 2 oc (n £ ). (9) 

Theorem 4 (1) The spectrum of the operator H £ consists of the absolutely continuous component 
[Ao,oo) (if n e has at least one infinite channel) and, possibly, a discrete set of positive eigenvalues 
X = fij i£ with the only possible limiting point at infinity. The multiplicity of the a. c. spectrum changes 
at points X = X n , and is equal to m(n + 1) 077 the interval (A„, A n +i). 

If VL e is a spider domain, then the eigenvalues p,j <e — fij do not depend on e. 

(2) The operator admits a meromorphic extension from the upper half plane Imk > into 
lower half plane Imk < with the branch points at k — zt^/X^ of the second order and the real poles 
at k = ±^//ij )£ and, perhaps, at some of the branch points (see the remark below). The resolvent 



flj) has a pole at k — ±yA~^ if and only if the homogeneous problem f3J) with A = A has a nontrivial 
solution u such that 

u = ajip n (y/e) + (e~ 7 *), x£Cj_ £ , t ~ > oo, 1 < j < m. (10) 

(3) If f G L 2 om (Q e ), and k — \f\ is real and is not a pole or a branch point of the operator (0|), 
and A > Ao, then the problem |[3J) ; is uniquely solvable and the outgoing solution u can be found 
as the Lf oc (fl £ ) limit 

u = R x+io f. (11) 

(4) There exist exactly m{n + 1) different scattering solutions for the values of A G (A„, \ n +i) 
such that k = \/A is not a pole of the operator and the scattering solution is defined uniquely 
after the incident wave is chosen. 

Remark. The pole of R\ at a branch point A = A„ is defined as the pole of this operator 
function considered as a function of z = \/X — X n . 

Let us describe the asymptotic behavior of scattering solutions ^S> = as e — + 0, A G (Ao, Ai). 
We shall consider here only the first zone of the absolutely continuous spectrum, but one can find 
the asymptotics of in [5T] for any A > Ao- Note that an arbitrary solution u of equation (fT|) in 
a channel Cj. £ can be represented as a series with respect to the orthogonal basis {ip n (y/e)} of the 
eigenfunctions of the Laplacian in the cross-section of Cj >e . Thus, it can be represented as a linear 
combination of the travelling waves 



e ±l - nt (p n (y/e), Ag(A„,A„+i), 

and terms which grow or decay exponentially along the axis of Cj, e . The main term of small e 

asymptotics of scattering solutions contains only travelling waves, i.e. functions \E , p E ' ) in each channel 
Cj, e have the following form when A G (Ao, Ai) : 

* = ^ = (c^e-*^' + /3 pje l ^^>o(y/e) + x G Cj, e> (12) 

where 

-ld(t) 

\r e p j\<Ce * , 7 > 0, and d(t) = mxn(t,lj — t). 
The constants ctp,j, (3 p .j and functions r p j depend on A and e. Formula (|12p can be written as follows 

* = w (y/ E )+r;, Ir^KCe- 2 ^ 1 , (13) 

where the function <; = ?(£) can be considered as a function on the limiting graph T which is equal 
to ?j(t), < t < lj, on the edge Tj and satisfies the following equation: 

d 2 

(e 2 ^2 +A-A ) ? = 0. (14) 

In order to complete the description of the main term of the asymptotic expansion (|12[) . we 
need to provide the choice of constants in the representation of ^ as a linear combinations of the 
exponents. We specify <; by imposing conditions at infinity and gluing conditions (GC) at each 
vertex v of the graph T. Let V = {v} be the set of vertices v of the limiting graph F. These vertices 
correspond to the free ends of the channels and the junctions in f2 e . 

The conditions at infinity concern only the infinite channels Cj jS , j < m. They indicate that the 
incident wave comes through the channel C p , e . They have the form: 

Pp,3= S P,j' ( 15 ) 
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The GC at vertices v of the graph T are universal for all incident waves and depend on A. We 
split the set V of vertices v of the graph in two subsets V = V\ U V2, where the vertices from the set 
V\ have degree 1 and correspond to the free ends of the channels, and the vertices from the set V2 
have degree at least two and correspond to the junctions Jj. e . We keep the same BC at v € Vi as 
at the free end of the corresponding channel of il £ : 

B e ( = atveVi. (16) 

In order to state the GC at a vertex v £ V2, we choose the parametrization on T in such a way that 
t = at v for all edges adjacent to this particular vertex. The origin (t = 0) on all the other edges 
can be chosen at any of the end points of the edge. Let d = d(v) > 2 be the order (the number of 
adjacent edges) of the vertex v € V 2 - For any function q on T, we form a vector ?W = $w(t) with 
d(v) components equal to the restrictions of <r on the edges of T adjacent to v. We shall need this 
vector only for small values of t > 0. Consider auxiliary scattering problems for the spider domain 
n„ e . The domain is formed by the individual junction which corresponds to the vertex v, and all 
channels with an end at this junction, where the channels are extended to infinity if they have a 
finite length. We enumerate the channels of fl Vi£ according to the order of the components of the 
vector We denote by T v the limiting graph defined by £l v , e . Definitions [IJ [2] and Theorem [4] 
remain valid for the domain Cl V)£ . In particular, one can define the scattering matrix T = 71, (A) for 
the problem {I} in the domain f2« jE . Let I v be the unit matrix of the same size as the size of the 
matrix D V (X). The GC at the vertex v € Vi has the form 

ie[I v +T v (X)}j t &Ht)-^X^Y [I v -T v (X)]&\t)=0 7 t = 0, veV 2 . (17) 

One has to keep in mind that the self-similarity of the spider domain fly e implies that T v = T V (X) 
does not depend on e. 

Definition 5 A family of subsets 1(e) of a bounded closed interval I C R 1 will be called thin if, for 
any S > 0, there exist constants /3 > and cy, independent of 5 and e, and ci = ci($), such that 1(e) 
can be covered by c\ intervals of length S together with c 2 e' 1 intervals of length C2e~@' e . Note that 
\l(e)\ -► as e -► 0. 

Theorem 6 For any bounded closed interval I C (Ao,Ai), there exists 7 = 7(/) > and a thin 
family of sets 1(e) such that the asymptotic expansion A13\) holds on all (finite and infinite) channels 
Cj <£ uniformly in A G I \ 1(e) and x in any bounded region of R d . The function q in US]) is a vector 
function on the limiting graph which satisfies the equation |j^[ ), conditions H15\) at infinity, BC I116\) , 
and the GC (77]). 

Remarks. 1) For spider domains, the estimate of the remainder is uniform for all x £ R d . 

2) The asymptotics stated in Theorem [6] is valid only outside of a thin set 1(e) since the poles of 
resolvent ([8]) may run over the interval (Ao, Ai) as e — * 0, and the scattering solution may not exist 
when A is a pole of the resolvent. These poles do not depend on e for spider domains, and the set 
1(e) is e-independent in this case. 

Consider a spider domain Cl v<£ and scattering solutions W = in fl v . £ when A belongs to a 
small neighborhood of Ao, i.e. 



* = = (Spje-*^-* + i PJ (A)e I ^ t )^ (2//£) + r s pd , x e C jt£ , |A - A | < 6, (18) 

where • decays exponentially as t — > 00 and t — t(x) is the coordinate of the point x. We define 
these solutions for all complex A in the circle | A — Ao | < 5 by the asymptotic expansion (| 18[) when 
ImA > 0, A 7^ A , and by extending them analytically for other values of A in the circle. 
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Lemma 7 Let fi„ iE be a spider domain. Then there exist 5 and 7 > such that 

1) for each p and |A — Ao| < 5, A 7^ Ao, the scattering solution exists and is unique, 

2) the scattering coefficients t p j{X) are analytic in \/\ — Ao when |A — Ao| < <5, 

3) the following estimate is valid for the remainder 

KJ < , A _ Ao | e ~^' ^ >0 ' xeCj, e . 

This statement can be extracted from the text of our paper [21] . Since it was not stated explicitly, 
we shall derive it from the theorems above. In fact, since the spider domain il v>e is self-similar, it 

is enough to prove this lemma when e — 1. We omit index e in Q, E , C p , £ , ty p when the problem in 
il e is considered with e = 1. Let a p (x) be a C°°- function on £1 which is equal to zero outside of 
the channel C p and equal to one on C p when t > 1. We look for the solution & p of the scattering 
problem in the form 



*p = <5 PJ e- lv ^^Vo(y)a P (x) + u p . 3A > 0, A ^ A , 
Then u p is the outgoing solution of the problem 

(-A -X)u = f, xeVL; u = Q on dfl, 

where 

/ = -5 pj [2V(e 4 ^ = ^Vo(y))Va p ( 2 :) + e' 1 ^™ Mv) ^a p (x)} € L 2 com (n). 

From Theorem2]it follows that there exists S > such that u p exists and is unique when | A — Ao | < 
3A > 0, A ^ Ao, and u p — R\f when 3A > 0, and u p can be extended analytically to the lower 
half-plane if R\ is understood as in @. The function u, p = R\f may have a pole at A = A . In 
particular, on the cross-sections t = 1 of the infinite channels Cj, the function u p is analytic in 
VA — Ao (when | A — Ao | < 6) with a possible pole at A = A . We note that / = on Cj fl {t > 1}. 
We represent u p there as a series with respect to the basis {ip s (y)}. This leads to (fT8|) and justifies 
all the statements of the lemma if we take into account the following two facts: 1) the resolvent (J9j) 
can not have a singularity at A = Ao of order higher than 1/|A — Aq|, since the norm of the resolvent 

at any point A does not exceed the inverse distance from A to the spectrum, 2) the scattering 
coefficients can not have a singularity at A = Ao due to Theorem [3J 

The proof of Lemma [7] is complete. 

4 Spectrum of the problem on the limiting graph. 

Let us write the inhomogeneous problem on the limiting graph T which corresponds to the scattering 
problem (fT4f , (fT5|) , (fT6| , (fT7|) . We shall always assume that the function / in the right-hand side in 
the equation below has compact support. Then the corresponding inhomogeneous problem has the 
form ^ 

(£ 2 ^ + A-AoK = / onT, (19) 



B e s = atveVi, ie[I v + T v (A)] -? M (0) —\/X—Xq [I v - T v (A)]?< B > (0) = at v e V 2 , (20) 

at 



$ = foe 1 ' , t » 1, on the infinite edges Tj 1 < j < m. (21) 

This problem is relevant to the original problem in f2 e only while A < Ai, since more than one 
mode in each channel survives as e — + when A > Ai. The latter leads to a more complicated 
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problem on the limiting graph (see 21 ). We are going to use the problem (fTT))) - (|2"Tj) to study the 
spectrum of the operator H £ when 



A = A + £>, |/x|<c. (22) 

As we shall see later, if f2 £ has a channel of finite length, then the operator H £ has a sequence of 
eigenvalues which are at a distance of order 0(e 2 ) from the threshold A = Ao. For example, if fl £ 
is a finite cylinder with the Dirichlet boundary condition (see the introduction) these eigenvalues 
have the form Ao + e 2 m 2 /l 2 , m > 1. Thus, assumption (I22| allows one to study any finite number 
of eigenvalues near A = Ao. 

Let us make a substitution A = Ao + e 2 fi in (fT9"l) - (f2"Tj) . Condition ([2H1) may degenerate at fi = 0, 
and one needs to understand this condition at /i = as the limit when [i — > after an appropriate 
normalization which will be discussed later. 

Lemma 8 . There is an orthogonal projection P = P{v) in R d ( v ' such that the problem US^ - Hil]) 
under condition H2S\) can be written in the form 

(^+MK = ^ 2 / onT- (23) 

B e s = atveVu Ps {v H0)+O{e)^ v) {0) = 0, P^ ^ {v) (0)+O(e)^ v) (0) = at v e V 2 , 

dt dt 

(24) 

q = (ije t ^ it , t»l, on the infinite edges T 3 l<j<m, (25) 

where d(v) is the order (the number of adjacent edges) of the vertex v and P 1 - is the orthogonal 
complement to P. 

Remarks. 1) The GC ([2"3")) at v € V2 looks particularly simple in the eigenbasis of the operators 
P, P^-. If e = 0, then it is the Dirichlet/Neumann GC, i.e after appropriate orthogonal transformation 

£l (0) = • • ■ = a (0) = 0, (0) = ■ • • = ^ (0) = 0, k = RankP. 

2) We consider fi as being a spectral parameter of the problem (|23|) - ([25l) . but one needs to keep in 
mind that the terms 0{e) in condition (|23[) depend on /1. 

Proof. Let us recall that the matrix T V (X) is analytic in y/X — X due to Lemma [7] Theorem |3] 
implies the existence of the orthogonal matrix C v such that D v :— C~ 1 T(Xq)C v is a diagonal matrix 
with elements v s — ±1, 1 < s < d(v), on the diagonal. In fact, from Theorcm[3]it follows that, for any 
A 6 [Ao, Ai], one can reduce T(A) to a diagonal form with diagonal elements v s — v s {X) where \ v s \ = 1. 
Additionally, one can easily show that the matrix T{Xq) is real-valued, and therefore, v s = ±1 when 
A = Ao. The statement of the lemma follows immediately from here with P = 4(7 — D V )C~ . 

Consider the Green function G e (pf, 70, fi), 7,70 € T, 70 ^ V, of the problem (|2"51) - ([2"5")) which is 
the solution of the problem with £~ 2 / replaced by 5 7o (7). Here 5 7o (7) is the delta function on T 
supported at the point 70 which belongs to one of the edges of V (70 is not a vertex). The Green 
function G (7,7o,m) is the solution of p3 l) -([2g )l with e = in ([24)). 

Let us denote by Pq the closure in L 2 (T) of the operator —j^z defined on smooth functions 
satisfying (|24p with e = 0. Note that the conditions (f2"4")l with s = do not depend on fi. Hence, P 
is a self-adjoint operator whose spectrum consists of an absolutely continuous component {// > 0} 
(if r has at least one unbounded edge) and a discrete set {fij} of non-negative eigenvalues. Let us 
denote by B c the disk < c of the complex /i-plane. 
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Lemma 9 For any c > there exist Sq > such that 

1) the eigenvalues {/i s (e)} of the problem i23\) -l[2~5 \) in the disk B c of the complex /i-plane are 
located in C\£ -neighborhoods of the points {fJ-j}, C\ — Ci(c), and each such neighborhood contains 
p eigenvalues /i s (e) with multiplicity taken into account where p is the multiplicity of the eigenvalue 

2) the Green function G e exists and is unique when \x € P c \{/x s (e)} and has the form 

r , * s 9(l, To, H,s) 

G e (7,7o,M) = — T7 s — , 

h(fi,e) 

where g is a continuous function 0/76 T m , 70 6 r„, [i € B c , e € [0,£q\, functions g and h are 
analytic in ^fjl and e, and d has zeros in B c only at points fi — /i s (e)- Here T m , T n are arbitrary 
edges of T. 

Proof. We denote by t = t(j) (i = t (j )) the value of the parameter on the edge of T which 
corresponds to the point 7 G T m (70 6 r„, respectively). We look for the Green function in the form 

G e =S m , n — — +a m ,»e^ t + 6 m .» e -^ t , 7 € T m , 70 6 T n , (26) 

where 5 m ^ n is the Kronecker symbol, and the functions a m ni b m n depend on £ , Mi £ - Obviously, (f2"5]) 
with e~ 2 f replaced by 5 70 (7) holds. Let us fix the edge T n which contains 70. We substitute (f2"6"]) 
into (|24p . (|23|) and get 2N equations for 2N unknowns a mj „, &,„,„, 1 < m < N, n is fixed. The 
matrix M of this system depends analytically on ^/JI, fi € B c and e 6 [0,£o]- The right-hand side 
has the form c\e % ^ to + C2e~ l ^ to , where the vectors ci, C2 depend analytically on ^JJi and e. This 
implies all the statements of the lemma if we take into account that the determinant of M with 
e = has zeroes at eigenvalues of the operator Po- The proof is complete. 

In order to justify the resolvent convergence of the operator H e as e — ► and obtain the asymp- 
totic behavior of the eigenvalues of the problem ( Q} near Ao we need to represent the Green function 
G e of the problem on the graph in a special form. We fix points 7* strictly inside of the edges Tj. 
These points split T into graphs F™* which consist of one vertex v and parts of adjacent edges up 
to corresponding points 7,. If T v is the limiting graph which corresponds to the spider domain £l v , £ , 
then r° u< is obtained from T v by cutting its edges at points 

When e > is small enough, equation (|23| on T v has d(v) linearly independent solutions satisfy- 
ing the condition from (|24|) which corresponds to the chosen vertex v. This is obvious if e = (when 
the components of the vector C" 1 ^^ satisfy either the Dirichlet or the Neumann conditions at v). 
Therefore it is also true for small e > 0. We denote this solution space by S v . Let us fix a specific 
basis {^,11(7, e)}, 1 < P < d(v), in S v . It is defined as follows. Let us change the numeration of 
the edges of T (if needed) in such a way that the first d(v) edges are adjacent to v. We also choose 
the parametrization on these edges in such a way that t = corresponds to v. Then 

ip P ,v = s^e-^ + t M (\y^\ 7 e iy 

Here t — t(j), X = X + e 2 fi, t p j = t^j are the scattering coefficients for the spider domain Q Vi£ . 
Obviously, i/j pv satisfies conditions (|24|) , and formula (|18|) can be written as 

®p}v = ^(7^ ; e)vo(y/e) +t e PiV x € C* i:£ , 7 = 7(2;). (27) 

where 7(2;) € Tj is defined by the cross-section of the channel Cj i6 through the point x. 

We shall choose one of the points jj in a special way. Namely, if 70 G T n then we chose 7„ =70- 
Then G £ belongs to the solution space S v , and from Lemma [9] we get 
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Lemma 10 The Green function G e can be represented on each part r° of the graph T in the form 

G e {l,lo,v) = tj- x ^2 Op,«^p,«(7>M> E )> M e -B c , e € [0,e ], (28) 
i<p<«i(«) 

where the function h is defined in Lemma\Q and a p>v — a Pi ^(7o, /Lt, e) are continuous functions which 
are analytic in e and yfjl 



5 Resolvent convergence of the operator H E . 

We are going to study the asymptotic behavior of the resolvent R\_ E — (H e — A) -1 of the operator 
H e when ([22]) holds and e — ► 0. When /j, is complex, the resolvent R\ is understood in the sense 
of analytic continuation described in Theorem [J] In fact, we shall study R\f only inside of the 
channels Cj j£ and under the assumption that the support of / belongs to a bounded region inside of 
the channels. We fix finite segments C Tj of the edges of the graph large enough to contain the 
points 7j. Let V = UP^-. We denote by C' e — UCj e the union of the finite parts C'^ e of the channels 
which shrink to P^- as e — > 0. We shall identify functions from L 2 (C' £ ) with functions from L 2 (Q, £ ) 
equal to zero outside C' E . We also omit the restriction operator when functions on fi e are considered 
only on C' s . 

K/S L 2 {C' £ ), denote 

^ £ ) = 1| / . , ^ ( /!(n )> = / fMv/e)dy/\\My/e)\\ L ,, 7er. 

We shall use the notation G e for the integral operator 

GeCh) = J G £ (7,7o,M)C(7o)rf7o, 7 G r - 

Theorem 11 Let /ioM. TTien /or any disfc -B c , i/iere exist £o = £o( c ) and a constant C < oo 
smc/i i/iat t/ie function 

Rx„ e f = (H £ - A)- 1 /, feL 2 (C' £ ), 

is analytic in ^JJl when ji € B c \0(e), where 0(e) is Ce -neighborhood of the set {fJ-j}, and has the 
form 

\\R x , e f-ip {y/e)Gj^,e)\\ L , { c.) < Ce\\f\\ L , {c ,y 

Remarks. 1) The points /ij were introduced above as eigenvalues of the problem (|23|) - ([25|) 
on the graph with e — 0. The GC in this case are the Dirichlet and Neumann conditions for the 
components of the vector C~ 1 s t - V ' . Obviously, these points are also eigenvalues of the operator H £ 
with the junctions of fl e replaced by the same Dirichlet/Neumann conditions on the edges of the 
channels adjacent to the junctions. 

2) The resolvent convergence stated in the theorem implies the convergence, as e — > 0, of eigen- 
values of operator H £ to {nj}- We could not guarantee the fact that the eigenvalues of the problem 
on the graph are real (see Lemma [9]). Of course, they are real for operator H £ . 

Proof. We construct an approximation K\ £ to the resolvent R\ £ = (H e — A) -1 for / e L 2 (C' £ ). 
We represent L 2 (C' £ ) as the orthogonal sum 

L 2 (C>) = L 2 (C>)+L 2 (C>), 

where functions from Lq(C' £ ) have the form h(j)ipo(y/£), 7 G V, and functions from L\(C' e ) on each 
cross-section of the channels are orthogonal to (po(y/e). Here and below the point 7 = 7(2:) € T is 
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defined by the cross-section of the channel through x. We put 70 = 7(0:0), i- e - 7o is the point on the 
graph defined by the cross-section of the channel through xq. 
Consider the operator 

K x . £ :L 2 (C' s )^L 2 (n e ) 
with kernel K X}£ (x, Xq) defined as follows: 

K\, e (x,x ) = ^2 r? 1 \ ^2 a P,v^ p e l(x,x ), X £ C £ , X e fi e . 

veV ^ £, ^> l<p<d(v) 

Here h(/i,e) and a P) „ — a p , v (70) Mj £ ) are functions defined in (|28[) . and Vpp^ are defined by the 
scattering solutions ^fpj, of the problem in the spider domain f2„ jE in the following way. Let e be 
the part of the spider domain Q„ e which consists of the junction and parts of the adjacent channels 
up to the cylinders C'j £ . Let Qf v e (fl^ £ ) be a bigger domain which contains additionally the parts of 
the cylinders C'j £ up to the cross-sections which correspond to points (the whole cylinders C'j E , 

respectively). We put ^p,l = ^p,l in fiJJ.e- We s P nt the scattering solutions in the cylinders 
C'j e into the sum of two terms. The first term contains the main modes ipo(y / 'e)e ±l ^ lt , and the 
second one is orthogonal to <po(y/e) in each cross-section. We multiply the first term by the function 
6 v (x, xq) equal to one on Q' v E and equal to zero everywhere else on Q £ . We multiply the second term 
by an infinitely smooth cut-off function r] v (x) equal to one on fl® e and equal to zero on fl e outside 
fij e . In other terms, 

$l c } v (x,x ) = e v (x,x )^l + {r] v {x) - e v (x 1 x )y pv (29) 

where r z V j is defined in ([27]) . 

Recall that the representation L = ur™* depends on the choice of points 7 S € T' s C T s . All these 
points are fixed arbitrarily except one: if 70 € Tj then 7^ is chosen to be equal to 70. This is the 
reason why 9 V depends on xq and r\ v is xo-independent. 

We look for the parametrix (almost resolvent) of the operator H e , when (j22[) holds and / S 
L 2 (C^), in the form 

F Xt£ = K XtS P + R X ,ePl, 

where Pq, Pi are projections on Lq(C' £ ) and Lf(C' £ ), respectively. It is not difficult to show that 
||i? A , e Pi|| = 0{e 2 ) and H £ F x , e =I + F x , e , where ||F A ,e|| = 0(e). This implies that R X<E = K Xt£ P Q + 
O(e). The latter, together with Lemma [3 justifies Theorem [TT1 

6 The GC at A near the threshold Ao- 

Theorem [TT] and the remarks following the theorem indicate that the GC at each vertex when 
A — Ao = 0(e 2 ) is the Dirichlet/Neumann condition, i.e. the junctions of £l £ can be replaced by k(v) 
Dirichlet and d(v) — k(v) Neumann conditions at the edges of the channels adjacent to the junctions 
(after an appropriate orthogonal transformation). We are going to specify the choice between the 
Dirichlet and Neumann conditions. First, we would like to make four important 
Remarks. 1) Classical Kirchhoff's GC corresponds to k = d — 1. 

2) For any domain fl £ under consideration, if A — Xq = 0(e 2 ) and the Neumann boundary 
condition is imposed on dfl e (Ao = in this case) then the GC on the limiting graph is Kirchhoff's 
condition (see section 2). 

3) It was proven in [21] that if A — Ao = 0(e 2 ) and the boundary condition on dVL £ is different 
from the Neumann condition, then the GC on the limiting graph is the Dirichlet condition (k = d) 
for generic domains Q £ . An example at the end of the next section illustrates this fact. 



12 



4) The theorem below states that Kirchhoff's GC condition on the limiting graph appears in the 
case of arbitrary boundary conditions on dQ e , if the operator H E has a ground state at A = Ao- The 
ground state at A = Ao = exists for an arbitrary domain f2 e , if the Neumann boundary condition 
is imposed on dfl E . The ground state at A = Ao does not exist for generic domains fi E in the case of 
other boundary conditions (see [21]). 

Note that the GC is determined by the scattering matrix in the spider domain Cl v , e , and this 
matrix does not depend on e. Thus, when the GC is studied, it is enough to consider a spider 
e-independent domain Cl = fi^i. We shall omit the indices v and e in H e , Cj. e when e = 1. 

Definition 12 A ground state of the operator H in a spider domain Q at A = Ao is the function 
"00 = i>o (%)> which is bounded, strictly positive inside ft, satisfies the equation (—A — \o)ipo — in 
ft, and the boundary condition on dft, and has the following asymptotic behavior at infinity 

V> (a) = ¥>o (l/) [pj +0(1)], x G C 3 , \x\ -> +oo, (30) 

where pj > and ipo is the ground state of the operator in the cross-sections of the channels. 

Let us stress that we assume the strict positivity of pj. 
Let's consider the parabolic problem in a spider domain fl e , 

du E y 

= Au e , u e (0,x) =<p(7)Vo(-), (r,x) \on e = 0, (31) 

where 7 = 7(2;) G T is defined by the cross-section of the channel through the point x, function ip is 
continuous, compactly supported with a support outside of the junctions, and depends only on the 
longitudinal ("slow") variable t on each edge Fj C V. We shall denote the coordinate t on Cj and 
Tj by tj . Let u>' be a compact in the cross-section u> of the channels Cj . 

Theorem 13 . Let Q be a spider domain, the Dirichlet or Robin boundary condition be imposed 
at dd, and let the operator H have a ground state at A = Ao. Then asymptotically, as e — > 0, the 
solution of the parabolic problem \31\) in fl e has the following form 

— x °j y 

u e (T,x)=e ~iz~w e {T,x)cpo(-), 

e 

where the function w e converges uniformly in any region of the form < c _1 < r < c, tj(x) > 6 > 0, 
y G lu 1 to a function w(t, 7) on the limiting graph T which satisfies the relations 

dw d 2 w \ dw 

on I ; w is continuous at the vertex, /^PjT; — (0) = 0. (32) 



dr dt 2 ' ' ^ ' 3 at 
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Remarks. 1) Let's note that under the ground state condition, operator H e has no eigenvalues 
below Ao- Otherwise, the eigenfunction with the eigenvalue A m i n < Ao must be orthogonal to the 
ground state ipo (x) , and this contradicts the positivity of both functions. 

2) The eigenvalues below Ao can exist if H does not have the ground state at Ao- For instance, 
they definitely exist if one of the junctions is "wide enough" (in contrast to the O. Post condition 
[24]). The solution u e (t,x) in this case has asymptotics different from the one stated in Theorem 
ri3l In particular, if the function ip (see ([31]) ) is positive, then 

e 2 lnu £ (t, x) — > A min , r — > 00. 

What is more important, the total mass of the heat energy in this case is concentrated in an 
arbitrarily small, as e — * 0, neighborhood of the junctions. The limiting diffusion process on Y 
degenerates. 
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Proof. For simplicity, we shall assume that the Dirichlet boundary condition is imposed on 
dVl £ . Obviously, the function ipa is the ground state in the spider domain il £ . In particular, 

£ 2 A^>o + A ^o = 0; V'o (^) =</?o(|)[ft + o(l)], xeC w \x\ -> +00; ipo\en c = . 

(X\ A T 
—J e~~^ r w £ (t,x). Then 

^ = Aw e + ^V (lnipo(^)) -Vw e , 
w e (0,x) = <p(i)6(^), 6(^)=j+o(l), xeCj, \x\^+oo. (33) 

We look for bounded solutions u S) w £ of the parabolic problems. We do not need to impose 
boundary conditions on dfl £ on the function w £ since the boundedness of w £ implies that u £ = 
on dft £ . The parabolic problem (|33p has a unique bounded solution (without boundary conditions 
on dfl £ ) since V 2 ln^o(-) is growing near dH, £ . This growth of the coefficient in (|33|) does not 
allow the heat energy (or diffusion) to reach dfl £ . The fundamental solution q £ — q s (r, Xo,x) of 
the problem (|33|) exists, is unique, and J n q £ (t,xo,x) dx = 1. This fundamental solution is the 

transition density of the Markov diffusion process — (T^ , Yp e ' ) in Q £ with the generator 

W e = A+^(vinVo(~),v) . 

Let Ti = H £ , and q = q £ when e = 1. The coefficients of the operator Ti are singular at the 
boundary of the domain. However, the transition density q{r, xo,x) is not vanishing inside fl. To 
be more exact, the Doblin condition holds, i.e. for any compact lu' C uj, there exist S > such that 
for any channel Cj the following estimate holds 

q(r, x ,x) > 5 when T > r > 1, x = (t,y), x = (t, y Q ), y, y Q E J . 

The operator Ti = A + 2 (V In ^0(^)1 V) has a unique (up to normalization) invariant measure. 
This measure has the density tt (x) = tpl (x) . In fact, if H = A + (VA, V) then H* = V - (VA, V) - 
(AA), and one can easily check that TL*e A ^ = 0. If we put now A (x) = 21n^ (a;)j we get 
n*n(x) = 0. 

When t > So > and e — > the transversal component and the longitudinal component 
T { T e) of the diffusion process in Q £ are asymptotically independent. The transversal component 

oscillates very fast and has asymptotically (e — > 0) invariant measure -tp^ (— ^. The latter follows 
from the Doblin condition. The longitudinal component has a constant diffusion with the drift which 
is exponentially small (of order 0(e - ^), 7 > 0) outside any neighborhood of the junction. 

Under conditions above, one can apply (with minimal modifications) the fundamental averaging 
procedure by Freidlin-Wcntzel (see [12] ) which leads to the convergence (in law on each compact 
interval in r) of the distribution of the process x| £ ' to the distribution of the process on T with 
d 2 

the generator — — on the space of functions on T smooth outside of the vertex and satisfying the 
dt A 

appropriate GC. The GC are defined by the limiting invariant measure. This limiting measure on 
r is equal (up to a normalization) to pj on edges Tj. This leads to the GC (|3"2"|) of the generalized 
Kirchhoff form. The proof is complete. 

Theorem 14 . Let operator H in a spider domain £1 with the Dirichlet or Robin condition at dfl 
have a ground state at X = Ao, and let A = Ao + 0(e 2 ). Then the GC |J7| j has the generalized 
Kirchhoff form: ( is continuous at the vertex v and Ylj—x Pj^r (0) = 0. 

This statement follows immediately from Theorem Q2] since it is already established that the GC 
has the Dirichlct/Ncumann form. 
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7 Effective potential. 



As it was already mentioned earlier, the GC (jTTJ) is A-dependent. The following result allows one 
to reduce the original problem in Q £ to a Schrodinger equation on the limiting graph with arbitrary 
A-independent GC and a A-independent matrix potential. The potential depends on the choice of 
the GC. Only the lower part of the a.c. spectrum Ao < A < Ai will be considered. It is assumed 
below that s = 1, and the index e is omitted everywhere. 

Let T (A) , A € [Ao, Ai], be the scattering matrix for a spider domain Q, let A_jy < A_jv+i < 
■ • ■ < A_i < A be the eigenvalues of the discrete spectrum of H below the threshold Ao- The 
function T (A) has analytic extension into the complex plane with the cut along [0, oo). It has poles 
at A — \-n, A-i. Let m_jv, • • • , m_i be the corresponding residues (Hermitian d x d matrices). 
These residues contain complete information on the multiplicity of Xj, j — —N, ■ ■ ■ — 1 and on the 
exponential asymptotics of the eigenfunctions ipj (x) , \x\ — > +oo. 

Theorem [6] allows one to reduce the problem in f2 e to an equation for a function £ on the limiting 
graph T with appropriate GC at the vertex. Consider the vector ip := Q v > = C t > 0, whose 

components are the restrictions of C to the edges of L. Note that the GC were formulated through 
the vector ^ V K Now we would also like to treat the equation for the function £ on T as a linear 
system for the vector on the half axis t > 0. 

Theorem 15 . There exists an effective fast decreasing matrix d x d potential V (i) such that 
V (t) = V* (t) , and the problem 

-ip" + [V{t)-X I}ip = Xil), t>0, ^(0) = (34) 

has the same spectral data on the interval (— oo, Ai) as the original problem in Q. The latter means 
that the scattering matrix S(X) of the problem coincides with T (A) on the interval [Ao, Ai], and 
the poles and residues of S{X) and T (A) are equal. 

Remarks. 1) The potential is defined not uniquely. 

2) The Dirichlet condition ip (0) = can be replaced by any fixed GC, say the Kirchhoff one (of 
course, with the different effective potential). 

3) Different effective potentials appeared when explicitly solvable models were studied in our 
paper [20] . 

Proof. This statement is a simple corollary of the inverse spectral theory by Agranovich and 
Marchenko for 1-D matrix Schrodinger operators [Tj. One needs only to show that J 1 (A) can be 
extended to the semiaxis A > Ai in such a way that the extension will satisfy all the conditions 
required by the Agranovich-Marchenko theory. 

Example to the statements of Lemma [8] and Theorem 1141 Consider the Schrodinger 

,2 

operator H — — 4p + v (t) on the whole axis with a potential v (t) compactly supported on [—1,1]. 
This operator may serve as a simplified version of the operator (|34[) . The simplest explicitly solvable 
model from [20] also leads to the operator H. The GC at t = for this explicitly solvable model 
are determined by the limit, as e — > 0, of the solution of the equation H £ ip e — f , where H e = 
— Jj? + £~ 2 v (I) and / is compactly supported and vanishing in a neighbourhood of t = 0. The 
solution ip s is understood as Lf oc limit of (H e + i[i)~ l f 6 L 2 as /i — * +0. 

Of course, A = is the bottom of the a.c. spectrum for H . If operator H does not have negative 
eigenvalues, then the equation Hip = has a unique (up to a constant factor) positive solution 
ipo (t), which is not necessarily bounded. If this solution is linear outside [—1, 1], then the limiting 
GC are the Dirichlet ones. This case is generic. If this solution is constant on one of the semiaxis, 
then the GC are the Dirichlet/Neumann conditions. Finally, if ipo (t) = p± for ±t > 1, then we have 
the situation of Theorem [Til the ground state and the generalized Kirchhoff 's GC. 

One can get a nontrivial Kirchhoff 's condition even in the case when H has a negative spectrum. 
It is sufficient to assume that A = is the eigenvalue (but not the minimal one) of the Neumann 
spectral problem for H on [—1,1]. 
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